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ABSTRACT 

The thermal infrared spectra of rocks are potentially of considerable 
importance for remote sensing. This paper presents results of a theoretical 
investigation of the relation between spectral features in the 8-12 microm- 
eter region and rock type. Data on compositions of a suite of rocks and 
measurements of their spectral intensities in 8.2-10.9 and 9.4-12.1 
micrometer bands published by Vincent (1973) were subjected to various 
quantitative procedures. It was found that there was no consistent direct 
relationship between rock group names and between the Thomton-Tuttle 
(1960) Differentiation Index and the relative spectral intensities. This 
relationship is explicable on the basis of the change in average Si-0 bond 
length which is a function of the degree of polymerization of the SK >4 
tetrahedra of the silicate minerals in the igneous rocks. 


v 


A PHYSICAL BASIS FOR REMOTE ROCK MAPPING 
OF IGNEOUS ROCKS USING SPECTRAL VARIATIONS IN 
THERMAL INFRARED EMITTANCE 


INTRODUCTION 

The use of infrared spectroscopy for the remote characterization of planetary and terrestrial sur- 
faces has received some interest due to efforts in the investigation of these bodies from space. The 
interested reader is referred to works such as those by Karr (1975) and Smith (1977) for some 
examples of these studies. 

Much of the emphasis in this work has been on the reflectance properties of rocks in the near infra- 
red (up to 2.5 jam) because the energy available in this reflective region of the spectrum is much 
greater than that at longer wavelengths. However, a depression in the thermal infrared emittance 
spectra of rocks (sometimes called restrahlen)is related to the sympathetic vibration of Si-0 bonds. 
Shifts in the location of this feature have been ascribed to variations in rock-type. This concept has 
spurred researchers such as Logan, Hunt and Salisbury (1975) and Lyon and Green (1975) to inves- 
tigate the relationships of the position of the emittance depression with respect to rock composition. 

In order to explore the utility of using this spectral variation for remote discrimination of rock 
types, NASA launched the Surface Composition Mapping Radiometer (SCMR) aboard the Nimbus-5 
spacecraft in 1972. This instrument had two spectral bands centered at 8.8 and 10.7 /un; variation 
in the relative intensities of these bands was to be correlated with variation in rock composition. 
Unfortunately, the instrument did not produce a great deal of data. However, there was some indi- 
cation that these data could be used for rock discrimination (Lyon and Green, 1975; Walter et al., 
1976). 
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Previous stu.! ivs < •! tv 'dationship between variations in spectral intensity in the thermal infrared 
and rock type have been largely empirical Vincent (1972) attempted to correlate the variation with 
silica content, but Si0 2 abundance by itself is notoriously poor as an indicator of rock type. He 
later presented a broader, but still empirically-based analysis of the problem (Vincent, 1973). We 
attempt here to derive statistically acceptable correlations between spectral variations and the varia- 
tions in the chemistry (and normative mineralogy) of rocks and to develop a rudimentary pheno- 
menological model for the variation in terms of mineralogical variations in rocks. From this, it is 
possible to develop a crude relationship between spectral variations and rock-type which may be 
useful for rock-classification by remote sensing from space. 

The work presented here is based on laboratory measurements of igneous ruck spectra published by 
Vincent (1973) together with corresponding chemical analyses and calculations of normative miner- 
alogy. The twenty-five samples analyzed covered a fiarly broad distribution of major igneous rock 
types, but unfortunately, for many of these types there were no replicate samples. The spectral 
variation was characterized as the ratio (R 12 ) of spectral intensities in a band from 8.2-10.9 urn 
relative to a band from 9.4-12.1 ptm. Sometime ago, Lyon (1964) showed that this ratio increases 
as the rocks become more basic (less silicic). 

Using these data, reproduced in Table 1, an attempt has been made to determine the major cause of 
this spectral change and relate this to some parameter which can be used to characterize igneous 
silicate rocks. 

Since the results of the analyses discussed below are used to make inferences to a larger population, 
it would be well for the reader to keep in mine that the 25 rock samples are igneous rocks. Further 
the proportions of the various rocks in the suite are not by implication considered to be representa- 
tive of any specific location on the earth’s surface. Vincent (1973) originally read this suite because 
of the wide range of textures and compositions represented. Therefore, the relationships developed 
within this paper are not necessarily applicable to metamorphic or sedimentary rocks. Proof of the 
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Tablt 1 

>ata Ua«d in This Work (from Vincont, 1 973) 
a) R«, > and Chomical Data (wt. %) 
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Table 1 (Continued) 

Normative Mineral Data (wt. %) (Cont’d) 
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Table 1 (Continued) 

b) Ri, > and Normativa Mineral Data (wt. %) 
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utility of procedures and conjectures described below for identifying such rocks as well as for 
specific spatial arrangements and proportions of rocks awaits further experimentation. 

RELATIONSHIP OF R| 2 TO ROCK NAME 

Before proceeding, it would be wise to determine whether or to what degree Rj 2 can be used to 
characterize, or act as a surrogate measure of, the relatively subjective classification given to hand 
specimens by a group of collectors. To test this hypothesis, the names bestowed upon the 25 igne- 
ous rocks by the geologist who collected them, were grouped as per Table 2. The analysis was per- 
formed by testing to see if there are significant differences in the Rj^ means of the groups. Accep- 
ting an hypothesis that claims there are significant differences in the means of the groups implies 
that the R^ value might be used to differentiate rocks which were largely identified in hand 
specimen. The procedure used to test hypotheses of this nature is called analysis of variance (Fisher, 
1970). Anorthite, diabase and peridotite were not used in this analysis because they represent one- 
member groups. The results of the analysis of variance (Table 3) were obtained using the MINITAB 2 
statistical package (Ryan et al., 1976). The F-Ratio is not significant at any conventional signifi- 
cance level and thus there is no evidence to indicate that these samples so grouped can be differen- 
tiated using Rj 2 * The large amount of overlap in the confidence limits illustrated in Figure 1 
graphically indicates the lack of significant difference between the means. 


Table 2 

Groupings of rock names for use in analysis of variance. 


Group 1 

Group 3 


Granite (A79) 
Granite (Al 19) 
Granite (A122) 
Rhyolite (SA-49) 

Granodiorite (A 1 1 7) 
Granodiorite (A 1 27) 

Group 4 

Andesite (E5SA) 
Andesite (2A3) 

Group 7 

Rhyolite 

Group 2 

Pyroxene Syenite (A83E) 
Trachyte 

Nepheline Syenite 

Dacite (E12A) 
Dacite (2A1) 
Dacite (2B5)W 

Group S 

Diorite (A129) 
Diorite (E47C) 

Gabbro (E80A) 
Basalt (1A8) 
Basalt (WI-1-104) 
Basalt (S95) 

Basalt (5P4) 

Basalt (1 A3) 
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Table 3 


Analysis of variance, test of significant difference among means of groups of samples. 


DUE TO 

DF 

SUM OF SQUARES 

MEAN SQUARE 

F* 

P(F > F*) 

GROUPINGS 

6 

0.008079 

0.001346 

2.19 

0.10 

ERROR 

16 

0.009832 

0.000615 



TOTAL 

22 

0.017911 






Group 


M * estimate of R, t mean 

*-—• » estimated 96% confidence limits 

Figure 1. Individual 95 Percent Confidence Interval for Group Means (Based on Pooled 
Standard Deviations). 
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The two related sets of data obtained by Vincent ( 1 972) were examined to determine the functional 


relationship (if any) between each of these data sets and the ratio R j ^ • The da'a sets, obtained from 
a suite of 25 igneous rock samples, were, 1) the weight percentage of 13 nutfor elemental oxides(here- 
after called “elemental data"); and 2) the weight percentage of 20 normative materials (hereafter 
called “mineralogical data"). The two data sets were not independently measured as the mineral- 
ogical data was derived from the elemental oxide data using a computer program devised by Rodgers 
et al., 1970. The values of the elemental data, the mineralogical data and R 12 for each rock sample 
have been reproduced in Table 1 . 

We have adopted the common statistical terminology designating the two compositional data sets 
as the independent or predictor variables and the ratio Rj 2 as the dependent or criterion variable. 
Through examination of each of the compositional data sets in turn, we attempted to establish a 
function which relates members of each of these sets to the corresponding criterion variable. These 
functions were to be parsimonious (Draper and Smith, 1966), i.e., have as few parameters as pos- 
sible, while being subject to the constraint of an adequate statistical fit and consistent with theore- 
tical considerations. 

The number of predictor variables must be reduced in order to 1 ) have sufficient degrees of freedom 
to perform a sensitive test of significance; 2) to avoid “artifacts" in the analysis such as overly opti- 
mistic estimates of diagnostic analysis statistics which can occur as the result of the inclusion of 
redundant or irrelevant variables (Draper and Smith, 1966); and 3) to reduce the influence of ratio 
correlations (Chayes, 1971 ) arising from having a predictor set whose values sum tc a constant (in 
this instance, each case sums to 100 percent). 
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The number of predictor variables was reduced in two steps. The first step used information which 
came strictly from the predictor data set; the objective being to reduce redundancy among the pre- 
dictor variables. The procedure involved extraction of components of the correlation matrix; pre- 
dictors which “load” strongly on the same component vary similarly across the sample space 
(Gorsuch, 1974). Because the method used for component extraction yields orthogonal compo- 
nents, variables loading highly on different components represented different patterns of variations 
and it is possible to reduce redundancy by selecting only one variable representing each pattern of 
variation. When there was a choice among more than one variable, theoretical considerations were 
used to guide the selection. The SPSS subprogram FACTOR (Kim, 1975) selecting options PA1 , 
and PA2 and a VARIMAX rotation was used as implemented on an IBM 370/3033 at the Pennsyl- 
vania State University. In the second step, the subset of predictor variables arising from the compo- 
nents analysis was input along with Rj^ values to an All-Possible-Subsets-Regression (Neter and 
Wisserman, 1974; Draper and Smith, 1966). This procedure evaluates regressions of the criterion 
upon all combinations of the predictor variables. There are ^ (^) regressions evaluated, where: p 

t 

is the number of predictors arising from step 1 , and (? ) * p'!/(p'-r)!r!. For example, in the case 
p' * 5, 31 regressions would be evaluated (5 with one predictor, 10 with two predictors, 10 with 
three predictors, 5 with four predict on, and one with all five predicton). The regressions are eval- 
uated by examining functions of the regression sum of squares (SSR), the error sum of squares 
(SSE), and plots of the residuals which are departures of the observed values of Rj 2 from their 
expected values. The program used to perform this analysis was the P9R program of the BMD sys- 
tem (Frane, 1979) using the Fumival and Wilson (1974) search procedure and R^j (defined as R 2 
adjusted for the number of degrees of freedom) ari-l Mallow’s Cp (Daniel and Wood, 1971 ) as the 
evaluation criteria. 
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E£syus 

Elemental Data 

For the elemental data, in step 1 , only phosphorous could be removed from consideration. The 
element had a greater than 80 percent (94.2) loading along with titanium (89.6) on component 1 . 
r Because silica represents the most dominant variable (having a mean of about 60% by weight), 

analyses were performed both with and without silica in the subset. This was done to determine 


what fine detail, if any, is being hidden by the presence of silica. 


For either subset (with or without silica) the R^j and C p criteria are both optimized when the 5 
elements Ca, K, Mg, Na, Fe +3 are the predictor variables. However, analysis of the residuals indicate 
that the peridotite value is an outlier and resulting in spuriously high coefficients of correlation 
(R*0.95) and determination (R 2 *0.90 and R^ D j*0.876) 


With peridotite removed from the analysis, the subset without silica achieves the lowest Cp value 
when Al, Mg and Fe +3 are in the equation. In the previously-cited analysis (with the peridotite 
Included), Al, Mg and Fe +3 provide the best regression if the number of independent variables used 
is restricted to three. Salient statistics from this regression equation are given in Table 4. 


Table 4 

Statistics for elemental data -- Si0 2 and P 2 Oj not included. 


Mallows’ Cp 0.S6 

Squared Multiple Correlation, R 2 0.685 

Multiple Correlation, R 0.828 

Adjusted Squared Mult. Corr., R 2 ^ 0.638 


Variable 

Name 

Regression 

Coefficient 

Standard 

Error 

Standard 

Coefficient 

T* 

P(T>|T*|) 

Contribution 
to R-Squared 

Intercept 

1.10993 

0.021564 

36.484 

51.47 

0.000 


Al 

0.00462 

0.001164 

0.500 

3.97 

0.001 

0.248 

Mg 

0.00844 

0.001689 

0.837 

4.99 

0.000 

0J92 

Fe + 3 

-0.01551 

0.006472 

-0.403 

-2.40 

0.026 

0.090 
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When silica is included in the analysis (no peridotite), it replaces A1 and Mg in the subset of inde- 
pendent variables forming the best regression equation (Table 5). 

The influence of Si in terms of explanatory power is equal to but opposite in “effect” to that of A1 
and Mg. Comparison of the results in Tables 4 and 5 indicates that the explanatory power (esti - 
mated by &adj) of the two models is roughly equal, while the coefficients of Fe 4 ^ in the two 
models are not significantly different from one another. Thus, the variation explained by Mg and A1 
is equal to that explained by Si. The inverse relationship between Al-Mg and Si can be seen in the 
opposite signs for the coefficients of Al-Mg (positive values) and Si (negative value). Consequently, 
the value of Rj^ decreases with increasing Si, or with decreasing A1 + Mg. Analysis of the residuals 
gives no evidence to reject the hypothesis that they are normally and independently distributed. 

Mineralogical Data 

Because there are 20 minerals in this predictor set with only 25 samples, the number of independent 
variables being considered must be reduced to permit formulation of a regression model. With the 


Table 5 

Statistics for elemental data — P 2 0 5 not included. 


Mallows’ C p -0.37 

Squared Multiple Correlation, R 2 0.646 

Multiple Correlation, R 0.804 

Adjusted Squared Mult. Corr., R 2 dj 0.613 


Variable 

Name 

Regression 

Coefficient 

Standard 

Error 

Standard 

Coefficient 

T* 

P (T>|T*1) 

Contribution 
to R- Squared 


1.39725 

0.367926 

45.929 

37.98 

0.000 


Si 

-0.00305 

0.000498 

0.963 

-6.13 


0.632 

Fe +3 

-0.01638 

0.006050 

-0.426 

-2.71 


0.124 







component - factor analysis of step 1, the number of predictors was reduced from 20 to 12. The 
variables selected are in bold face in the canonical form of the pattern matrix given in Table 6. 


Initial runs of the step 2 analysis, led to discarding nepheline syenite and peridotite as outliers be- 
cause of their large Mahalanobis D 2 values. Wollastonite had to be dropped from the predictor sub- 
set because it occurred in only the nepheline syenite sample. 

Two regression equations emerged from the analysis. The regression containing the independent 
variables forsterite, anorthite and water is statistically stable in that it emerges as the best (or 
nonsignificantly different from the best) subset regardless of the presence or absence from the sam- 
ple of nepheline syenite or peridotite. This regression also has the optimal Mallows’ Cp value of 
those models drawn from a subset of the predictor variables defined below. Summary statistics are 
given in Table 7. 


Table 6 

Factor pattern matrix for mineralogical data. 



FI 

F2 

F3 

F4 

F5 

F6 

F7 

CPWO 

0.909 







CPEN 

0.895 







CPFS 

0.879 







ILM 

0.920 







APA 

0.832 







MAG 

0.756 







QTZ 

-0.607 







ANOR 

0.531 







ORTHO 

OLFO 

-0.462 

0.983 



-0.648 



OLFA 


0.856 






ALB 

OPEN 


-0.644 

0.861 


0.514 



OPFS 

NEPH 



0.760 

0.945 




WOL 




0.983 




HEM 





-0.888 



CAL 






0.783 


COR 

WATER 






0.852 

0.883 
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Table 7 

Statistics from regression of water, forsterite and anorthite. 


Mallows’ Cp 2.56 

Squared Multiple Correlation, R 2 0.69776 

Multiple Correlation, R 0.83532 

Adjusted Squared Mult. Con., R 2 DJ 0.65004 


Variable 

Name 

Regression 

Coefficient 

Standard 

Error 

Standard 

Coefficient 

T* 

P(T>|T*|) 

Contribution 
to R-Squared 

Intercept 

1.17277 

0.00956622 


122.59 

0.000 


WAT 

-0.00452227 

0.00255464 

-0.230 

-1.77 

0.093 

0.049849 

OLO 

0.00232985 

0.00116707 

0.269 

2.00 

0.060 

0.063395 

AN 

0.00124317 

0.000256475 


4.85 

0.000 

0.373736 


The second predictor subset includes orthopyroxene, quartz, forsterite, clinopyroxene, oral 
(a summed combination orthoclase and albite) anorthite and water. It yielded a model containing 
oral, orthopyroxene and quartz (with the results for olivine leading us to an uncertain conclusion as 
to its significance). Statistics from this regression are given in Table 8. The order of “importance” 
of the variables, using the contribution to R 2 as the importance criterion, is quartz, alkali feldspar, 
pyroxene and olivine. 


Table 8 

Statistics from regression of oral (orthoclase+albite), orthopyroxene and quartz. 


Squared Multiple Correlation, R 2 0.609 

Multiple Correlation, R 0.780 

Adjusted Squared Mult. Corr., R 2 DJ 0.547 


Variable 

Name 

Regression 

Coefficient 

Standard 

Error 

Standard 

Coefficient 

T* 

P(T>iT*|) 

Contribution 
to R-Squared 

Intercept 

1.25502 

0.019532 

40.919 

64.25 

0.0 


ORAL 

-0.00090 

0.000347 

-0.479 

-2.61 

0.017 

0.140 

OPE 

-0.00213 

0.001476 

-0.247 

-1.44 

0.165 

0.043 

QTZ 

-0.00126 

0.000384 

-0.537 

-3.30 

0.004 

0.223 





INTERPRETATION 


The decreased emittance of silicate rocks in the region from 9 to 1 1 ftm is known to be due to the 
sympathetic vibration of the stretching of the silicon-oxygen bonds at these wavelengths. The 
frequency (at ca. 1000 cm -1 ) is directly related to the force by which the silicon and oxygen ions 
are bonded together; the greater the force, the higher the frequency of the sympathetic vibration 
and thus the shorter the wavelength at which one observes the decreased emittance. The oxygen 
of the Si-O couple is shared by another cation (Si-O-X) and the force of the Si-0 bond is affected 
by that of the O-X bond. Thus, the shape and location of the depression can change. 

Three types of effects on the Si-0 bond may be considered: substitution, disordering and polymer- 
ization. Substitution effects are the result of replacement of one cation by another in a single 
mineral species (e.g., Fe by Mg in orthopyroxene). Replacement of large cations by cations with 
smaller ionic radii causes the frequencies of the sympathetic vibrations of the Si-0 bond to decrease. 
This, in turn, results in an increase in the wavelength of Si-0 absorbtion (Fanner, 1974). [It is, 
perhaps, more appropriate to consider the relationship between the location of the emittance de- 
pression and the electronegativity of the secondary cation (Strens, 1974)] . Disordering of cations 
within potential crystal sites also has an effect. For example, random substitutions of aluminum for 
silicon in feldspar results in an emittance minimum at higher wavelengths than if the same feldspar 
is ordered (Kovach et al., 1975). Polymerization of silica tetrahedra can be considered an extreme 
condition of ordering. It increases from the case of isolated tetrahedra (e.g., olivine) through chain 
silicates (e.g., pyroxene) and culminates in those cases in which the silica tetrahedra are all bound 
together in a framework (e.g., quartz and feldspar). The effect is such that in the progression from 
isolated to framework silicates, the depression in spectral emissivity which occurs at approximately 
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10 jum decreases in wavelength. In order to relate this shift in the emittance depression to the com- 
position of minerals and rocks, it is first important to determine whether disordering, substitution 
or polymerization effects are the major cause of the shift. 

Disordering can be ruled out immediately ; this is an isochemical effect and thus cannot explain the 
variation of R 12 with chemistry as discussed in the previous section. 

Two of the major monomineralic substitutions in minerals are those between iron and magnesium 
(in pyroxenes and olivines) and aluminum and silicon (in feldspars). In the previous section, it was 
shown that, insofar as this suite of igneous rocks is concerned, there is a positive correlation between 
the amount of aluminum and magnesium in them and the values for R^. The centers of the two 
bands used to determine R^ lie at 9.55 and 10.75 Mm. The centers for the emittance depressions 
of these minerals (pyroxene, olivine and feldspar) lie between these band-center positions. Thus, 
for example, when ferrous iron substitutes for magnesium in a mineral, the wavelength of the emit- 
tance depression increases (Kovach et al., 1975) and the factor, Rj^ (the ratio of the intensity of 
the low-wavelength band to that of the high wavelength band) should increase. The reverse seems 
to be the case; Rj 2 is positively correlated with total rock magnesium content. Likewise, the sub- 
stitution of aluminum for silicon theoretically should reduce Rj 2 but we find that increased total 
aluminum content results in increased Rj^.* These relationships provide a strong Indication that 
solid substitution within mineral species is not the mechanism of the shift in Rj 2 . 

In chemical terms, the effect of polymerization on the position of the depression described above is 
consistent with the trend indicated in the previous section wherein Rj^ increases (reflecting an 

'This argument is based on the following assumptions which the authors consider to be warrented: 

1 . Substitutions of cations in mineral species mimics the variation of these cations in bulk rock chemistry; 

2. The shapes of the curves at the emittance depressions do not differ greatly from one rock type to another. 
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increase in the wavelength of the emittance depression) as aluminum and magnesium increase and 
silicon decreases. Thus, Rj 2 increases from granite to basalt to peridotite as the position of the 
emittance depression increases in wavelength. This is also consonant with the results of the statis- 
tical analysis which indicate that quartz, alkali feldspar, pyroxene and olivine, in that order, are 
related to variation in Rj^. Finally, this is consistent with the trend for the position of the emit- 
tance depression which increases from 9.2 to about 1 1 /an as we go from quartz to olivine. 

In crystallographic terms, polymerization results in changes in the Si-0 bond length and strength 
(Smith and Bailey, 1963) which, in turn, affects the R| 2 factor. Brown and Gibbs (1969) related 
bond length and strength to the average cation coordination of oxygen in silicate minerals. These 
authors determined a linear relationship between the average coordination number, CR and the 
mean Si-0 bond distance, Si-O: Si-O * 0.015 CM + 1 .579 (units given in Angstroms). Using this 
relationship, we can obtain values for the average Si-0 bond distance for the whole rock by weight- 
ing each mineral average by its abundance relative to other silicates in e?ch rock. The values for CN 
used for the major minerals were: quartz = 2.0; albite, orthoclase = 2.88 ; anorthite = 2.78; wolla- 
stonite, ortho- and clino-pyroxene = 3.33 ; olivine = 4.0 and nepheline = 2.83. These normative min- 
erals account for over 90% of every rock in this suite. The total weighted average oxygen coordina- 
tion number is plotted versus Rj 2 in Figure 2. Using all 25 samples for which Vincent ( 1 973) gave 
normative analyses, the correlation coefficient, r, was found to be 0.86. These and other statistics 2 


2 

The slope, m, and intercept, c, of the simple linear regressions reported in Tables 9 and 10 were calculated using a 
procedure which minimizes the sum of squared deviations from the x and y direction simultaneously. This avoids 
the common and inconect practice of first building a regression model which minimizes the y deviations for a 
given level of x and then using this model to predict x for a given level of y. The method of calculating m and c 
are from Cohen, 1980, and will be referred to as an invertible linear regression. 
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Figure 2. 



appear in Table 9. 3 The interpretation of these is that Si-0 bond distance, which is related to Si-O 
bond strength and polymerization of Si, is probably the m^jor cause of the shift in the position of 
the emittance depression. 


Table 9 

Invertible linear regression (using method of Cohen, 1980) between Rj ^ 
and average coordination number (CN): (Rj^ * roCN + c). 



m 

c 

r 

All 25 Samples 

0.177 

0.692 

0.86 

Without Sample 26 

0.127 

0.832 

0.68 

(Peridotite) 

Without Sample 24 

0.173 

0.703 

0.89 

(Anorthosite) 

Without Samples 24 and 26 

0.149 

0.768 

0.74 


3 The value of r in Tables 9 and 10 will vary depending on whether or not samples 24 (anorthosite) and 26 (perido- 
tite) are included in the analysis. The impact of these samples upon r can be seen by examining the following 
equations: i^*l -SSE/SSTO • SSR/SSTO, where r * v^andSSTO ” SSR + SSE. The position of these 
points (samples 24 and 26) relative to the rest of the points will determine which of the terms SSE or SSR (along 
with SSTO) will be increased. In the case of sample 26, because there are no other points in its range of x and y 
values, it will have a large impact on the regression line, causing the line to pass close to the point. Thus, while we 
are increasing SSTO, most of the increase is explained by the regression (SSR). This will produce an increase in 
P and hence in r. Alternatively, there are a number of x values in the same range as the x value of sample 24. 
Thus, the impact upon the regression of sample 24 is much diluted and so while there is an increase in SSTO, most 
of the increase is in the error term (SSE). This situation will produce a decrease in P and r. From a statistical 
viewpoint, the anorthosite is an outlier. Further appropriateness of the regression should not depend upon inclu- 
sion of the peridotite. The reader, however, is presented with the results of four regressions to allow consideration 
of all combinations for including samples 24 and 26. 
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CONCLUSIONS 


A theoretical basis for the empirical relationship between the position of the thermal infrared emit- 
tance depression exhibited by silicate igneous rocks has been established as being due to the charge 
in the average Si-0 bond strength and distance with increased average degree of polymerization c i 
SK >4 tetrahedra of the silicate mineral constituents. 


Interestingly, the trend of increased polymerization toward tectosilicates (quartz and feldspar) is 
the same trend observed in the crystal fractionation of igneous rocks. The evolution of rocks in a 
magmatic series is generally marked by increased concentrations of normative quartz and alkali feld- 
spars. On this basis, Thornton and Tuttle (1960) propos A the use of a quantitative measure of 
chemical differentiation which can be used for the classification of igneous rocks. The Differentia- 
tion Index (DI) is defined as the amount (in weight percent) of normative quartz + orthoclase + 
albite + nepheline + leucite + kalsilite. Thornton and Tuttle point out that, “no more than three of 
these normative minerals will appear in any given norm. Thus, the Differentiation Index is simply 
the sum of the percentages of the three normative minerals”. In applying the DI to chemical anal- 
yses for representative major igneous rocks, the authors obtained the following values: 


Intrusive Rock Type 

DI 

Extrusive Rock Type 

DI 

Alkali Granite 

93 

Alkali Rhyolite 

91 

Granite 

80 

Rhyolite 

88 

Granodiorite 

67 

Quartz Latite 

68 

Diorite 

48 

Andesite 

66 

Gabbro 

30 

Basalt 

35 

Olivine Gabbro 

27 

Olivine Diabase 

30 

Peridotite 

6 

Picrite 

12 


The Differentiation Indices for the rock suite used in this work are plotted against Rj 2 in Figure 3 
and statistics are given in Table 10. The DI’s for the major intrusive rock types are noted on the ordi- 
nate. (In view of the location of sample a, the peridotite, the use of a second-order equation would 
significantly improve the correlation. However, there is no theoretical basis to justify this step.) 
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Table 10 

Invertible linear regression (using method of Cohen, 1980) between differentiation index (DI) 

and Rj 2 -’ (R^ *mDI + c). 



m 

c 

r 

All 25 Samples 

-0.00172 

1.302 

-0.74 

Without Sample 26 

-0.00120 

1.270 

-0.68 

(Peridotite) 

Without Sample 24 

-0.00174 

1.309 

-0.72 

(Anorthosite) 

Without Samples 24 or 26 

-0.00116 

1.262 

-0.63 


The relationship on which Figure 3 is based appears to combine the fundamental causative principles 
behind the variation in the position of the emittance depression and those of igneous rock petro- 
genesis and classification. It may, therefore, be the most appropriate relationship for classification 
of rocks by remote sensing using this portion of the electromagnetic spectrum. 


One of the problems, however, in interpreting Vincent's data lies in the selection of the spectral 
bands used to determine Rj 2 ; the overlap between the bands undoubtedly weakens the sensitivity 
ofR u to the rock composition variable. Based on theoretical considerations and analyses of ther- 
mal infrared spectra of 134 rocks, Holmes and Nuesch (1978) proposed three improved spectral 
bands for rock discrimination in the thermal infrared region of the spectrum. These bands are : 

8.1 - 9.1 , 9.5 - 10.5 and 1 1 .0 - 12.0 pm. The use of these bands would probably improve upon 
the correlations obtained in this paper. 


It is interesting to note that the relationship between average Si-0 bond distance and the position 
of the emittance depression is independent of the texture of the rock. This is probably because the 
degree of polymerization is a statistical relationship based on the abundance of silicon relative to 
other cations; this relationship may well hold whether the material is crystalline or glassy. 
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While the evidence for the relationship between the position of the emittance minimum and the 
nverage Si-0 bond distance is rather compelling, the authors would feel more comfortable, espe- 
cially concerning the quantitative aspects of the relationships, if additional laboratory observations 
were available. We suggest that a significantly larger set of samples be analyzed. Indeed a great deal 
of the correlation we see rests on one sample, the peridotite, and there is a considerable hiatus from 
the next most basic rock (the gabbro at R J2 * 1.24) to the peridotite (R J2 * 1.372). Although 
rocks which fall between gabbroic and peridotitic compositions are rare, some samples should be 
found and analyzed to fill this void. 
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FIGURE CAPTIONS 


Figure 1 . Individual 95 Percent Confidence Interval for Group Means (Based on Pooled Standard 
Deviations). 

Figure 2. Average Oxygen Coordination Number (CN) Plotted Against Ratio of Intensities in Two 

Thermal Infrared Bands (Rj 2 ). Solid line all 25 samples; dashed line all 

samples except peridotite (a) and anorthosite (b). Si-0 distance, Si-O, calculated on the 
basis of CN is also given on the ordinate. 

Figure 3. Differentiation Index (DI) Plotted Against Ratio of Intensities in Two Thermal Infrared 
Bands (Rj 2 )• Solid all 25 samples; dashed line — all samples except perido- 

tite (a) and anorthosite (b). Intrusive rock types conesponding to DI values are indicated 


in the ordinate. 
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